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In this paper, we show that a Gaussian averaged neutrino oscillation model is equivalent to 
a neutrino decoherence model. Without loss of generality, the analysis is performed with two 
neutrino flavors. We also estimate the damping (or decoherence) parameter for atmospheric 
I . neutrinos and compare it to earlier obtained results. 
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I. INTRODUCTION 

^ ! 

CO , In this paper, we investigate two different models for transitions of neutrinos and how they are related to each 
' other. The two models are: neutrino oscillations and neutrino decoherence. Neutrino oscillation models have been 
I the standard description for neutrino transitions so far and still are. Neutrino decoherence models, on the other hand, 
^ • have recently been discussed by several authors [1-8] as an alternative description to neutrino oscillation models. 

[ There exist also other plausible descriptions of neutrino transitions such as neutrino decay models, which have been 
t — . suggested by Barger et al. [9-12]. However, these models will not be discussed in this paper. 

' The paper is organized as follows: In Sec. II, we go through the formalisms of neutrino oscillations and neutrino 
decoherence. In Sec. Ill, we give the condition for them to be equivalent with each other and we also try to estimate 
' ^ ' the decoherence term and in Sec. IV, we estimate the damping or decoherence parameter for atmospheric neutrinos. 
Finally, in Sec. V, we present the summary and also our conclusions. 

Ph! II. FORMALISM 

Oh: 

D . A. Neutrino oscillations 

The theory of neutrino oscillations is the far most plausible description of neutrino transitions. Neutrino oscillations 
between different neutrino flavors, Va and vp, occur with the well-known neutrino transition probabilities 
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where Sap is Kronecker's delta, L is the neutrino path length, E is the neutrino energy, n is the number of neutrino 
flavors, Am^jj = rn^ — m^, a, 6 = 1, 2, . . . , n is the (vacuum) mass squared difference between different neutrino mass 



*E-mail addresses: tohlssonOphysik. tu-muenchen. de , tonrniyStheophys . kth. se; Homepages: http://www.physik.tu-muenchen.de/ 
"tohlsson/ , http : //www. theophys .kth. se/~toinmy/ 



2 



eigenstates i^a and (or rather \i/a) and with masses nia and mt, respectively, and 



Ufil U^2 Ufj_3 ■ ■ ■ 
UtI Ur2 Ur3 ' ' ' 



V 



is the unitary n x n Maki-Nakagawa-Sakata (MNS) mixing matrix [13], which relates the flavor states l^a) {ct = 
e, /i, T, . . . ) and the mass eigenstates |z/a) (a = 1, 2, . . . , n). 

However, since in practice a neutrino wave is neither detected nor produced with sharp energy or with well-defined 
propagation length, wc have to average over the L / E dependence and other uncertainties in the detection and emission 
of the neutrino wave. 

We will here use the Gaussian average, which is defined by 



/oo 
P{x)f{x)dx, 
-OO 



(2) 



where 
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Here £ = (x) and a = ^J {{x — (x))^) arc> the expectation value and standard deviation, respectively. 

By taking the Gaussian average of Eq. (1) and using x = j^, we obtain the averaged transition probabilities from 
Va to Vjs as 

n n 
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The physical interpretations of the parameters £ and a are the following: 



The parameter £: The parameter i deals with the sensitivity of an experiment and is given by ^ = {L/E)/4. 
Note that we will here use {L/E) = {L)/{E), i.e., 



m' 



(4) 



This simplification holds if L and E are independent. 



• The parameter cr: The parameter cr is a so-called damping factor, which is responsible for the damping of 

the neutrino oscillation probabilities. 

A pessimistic upper bound for the damping parameter a is given by the uncertainty in a;, i.e., 
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where AL and AE are the uncertainties in the neutrino path length and the neutrino energy, respectively. A 
more optimistic upper bound would be 



< 



(6) 



For large values of cr, the dependence on the mass squared differences will be completely washed out, since 
1 - cos(2€ATO2(,)e-2'^'('^"«'')' 1 and sin (2€Ato2^) e-^'^'^^"*'")' when a ^ oo, and the Gaussian 
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averaged transition probabilities (Pa/3) will just be dependent on the MNS mixing matrix elements, the Uaa's, 



I.e., 



}v^JPc.0) = - 2 ^ Y\{UlJJpJJabU*p,). (7) 

a=l 6=1 

a<b 

Note that the imaginary part sum (the second sum in Eq. (3)) does not appear at all (in any form) in Eq. (7). 

Equation (7) corresponds to the classical limit. 

In the other limit, tr ^ 0, we will just regain Eq. (1) from Eq. (3) with I = x, i.e., 

lim(P„/3) = P^p. (8) 

(7 — ^0 

In the case of two neutrino flavors (n = 2), we can call these neutrino flavors Ve and u^, the unitary (orthogonal) 
2x2 MNS mixing matrix U is usually parameterized as 



TT - (TT \ - ( \ _ / COS 6* sin ^ \ 

U - {Uaa) - y jj^^ J - ^_sine cose J 



(9) 



where is the (vacuum) mixing angle. Note that U* = U in this case. Furthermore, Am^ = Am2x = —Am\2- 
Observe that for two neutrino flavors there cannot be any CV (or T) violation in the MNS mixing matrix U, since U 
is a real matrix in this case, which means that 

Pee = 1 - Pep = 1 - Ppe = PpM' (W) 

By taking the average of Eq. (10), we of course also have 

(Pee) = 1 - (Pe^) = 1 - (P^e) = (P^^)- (11) 

Thus, inserting a = e, (3 = ^, and n = 2 into Eq. (3), we obtain the Gaussian averaged two flavor neutrino oscillation 
formula as 

{PeM = \ Sin' 29 (l - e-2-'(^-')' cos (2Mm2)) , (12) 

where (. ~ {L)/{A{E)), i.e., 



{Pe,)m, {E}) = I sin^ 29 (l - e-^(--T cos ^) . 
In what follows, we will use (Pa/j) = Pa/3, {L) = L, and {E) = E for simplicity. 



(13) 



B. Neutrino decoherence 



Neutrino decoherence arises when we consider a neutrino system to be coupled to an environment (or a reservoir 
or a bath). In general, wo arc not allowed to use the Schrodiiigcr equation to describe system-environment typo 
interactions, since even if we start initially with a pure quantum mechanical state, the coupling to the environment 
will produce mixed quantum mechanical states. Thus, we are forced to use the Liouville equation. 

Assume that p = p{t) is the neutrino density matrix [21], which is Hcrmitian [p'^ = p) with positive eigenvalues and 
has constant trace equal to one (tr/o = 1). The Liouville equation for the neutrino density matrix is then 

p = -i[Hm,p], (14) 

where Hm is the free Hamiltonian for the neutrinos in the mass eigcnstatc basis. We will assume that the neutrino 
system is a two- level system, i.e., we have two neutrino flavors (n — 2), Vg and v^, as in the neutrino oscillation 
formalism. Thus, in this case, the neutrino density matrix /) is a 2 x 2 matrix and the free Hamiltonian Hm is also a 
2x2 matrix and is given by 

Hm = ^diag {ml, ml) or H'^ = Hm- ^(tr Hm)h = ^diag (-Am^, Am^), (15) 
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where I2 is the 2x2 identity matrix and again Am^ = — //tj^. 
of the Hamiltonian H^, since they just differ by a trace (of any of them). This trace will only give rise to a global 
phase that will not affect the transition probabilities anyway and is therefore irrelevant. 

Furthermore, note that we will not consider matter effects in this paper. In the case of neutrino transitions in 
matter, the Hamiltonian has a much more complicated structure. 

Solving Eq. (14) under the above assumptions leads to the ordinary neutrino oscillation transition probabilities 
in Eq. (1) with n = 2 as it should. However, the neutrinos could be influenced by so-called decoherence effects. 
Such effects are introduced by an extra term T>[p] in the Liouville equation for the neutrinos, i.e., the Markovian 
Liouville-Lindblad quantum mechanical master equation [2, 14, 15] 

p=-i[Hm,p]-V[p], (16) 

which allows transitions from pure quantum mechanical states to mixed quantum mechanical states. We will here use 
the Lindblad form for the decoherence term [14], viz., 

n 
a=l 

where the .Da's are Lindblad operators [22] arising from tracing or averaging away environment dynamics and must 
be such that '^^=1 ^a^a is a well-defined 2x2 matrix. Note that the extra term T^lp] is responsible for the fact that 
the quantum mechanical states can develop dissipation and irreversibility, and possible loss of quantum coherence i.e., 
(quantum) decoherence [2, 16]. The Da's are normally Hermitian {D^ = Da) if we require monotone time increase of 
the von Neumann entropy, S = — tr(plogp), and they usually also commute with the Hamiltonian {[Hm,Da] = 0) if 
we require conservation of the statistical average of the energy, i.e., -^tr (Hmp) = 0. 
Assuming that the D^s are Hermitian, Eq. (17) becomes 

n n 

1^[P] = E P^^ = E (P^l + ^Ip - '^DapDa) . (18) 

a— 1 a—1 

The operators which are all Hermitian, can therefore be expanded in the Pauli matrix basis as 

p = l(/2 + p.eT), (19) 

Hm = (20) 

Da = ^da-o-, a=l,2,...,n, (21) 



where I2 is again the 2x2 identity matrix, cr = (ui, C72, 0-3) is the Pauli matrix vector, and k = (0,0, —fc). Here 

IE 



k = ^^g- . In the Pauli matrix basis, we have after some tedious calculations 



P = ^P-*^' (22) 
[Hm,p] = ^(kxp)-o-, (23) 



1 " 

m = 2 E (Idal'P • - (da • P)(d„ • a)) , (24) 

a=l 

which means that Eq. (16) can be written in the form 

n 

p-cT = (kxp)-a-E (Idal'p • «T - (d„ • p)(d„ • a)) . (25) 

a=l 

Using well-known formulas from vector algebra, the sum in the above equation can be expressed as 

n 

E (da X (da X p)) • <r. 
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Thus, we have 

n 

p • cr = (k X p) • <7 - ^ (d„ X (d„ X p)) • cr (26) 



a=l 

or 



p = kx p-^da X (da X p) = k X p - ( X]|da|^ J p + ^(d^ • p)da, (27) 

a=l \a=l / a=l 

which is the Bloch-Lindblad equation. The first term after the second equality sign in the Bloch-Lindblad equation, 
k X p, is giving rise to neutrino oscillations, whereas the second and third terms are the decoherence terms. 

If [Hm, Da] = 0, then it follows that k x da = 0, i.e., k and da are parallel vectors. Thus, we can put da = dak, 
where k = k/|k|. Inserting this into Eq. (27), we obtain 

{Pi,P2,P3) = {kp2,-kpi,0) - (f{pi,P2,P3) + {0,0, (fps), (28) 

where cP = YZ=id'l^ 

Pi = kp2-(fpi, (29) 

P2 = -kpi- (fp2, (30) 

P3 = 0, (31) 

which can be written in matrix form as (p = Mp) 

/ k 0\ / pi\ 

(32) 





k 










•) 















This system of first order differential equations has the solution 

Pi{t)\ I e-'^^^coskt e-'^^^sinkt 0\ /pi(0) 

p{t)= \ P2{t) = -e-'^'*sinfct e-'^'*cosH P2{0) | = e^*p(0), (33) 

P3{t) V 1/ \pm 



p^{t) = pi(0)e-'^ *cosfct +P2(0)e-'' *sinfct, (34) 
P2(i) = -pi(0)e-'''*sinA;i + p2(0)e-'''*cosA;t, (35) 
P3{t) = P3{0). (36) 

Another, but equivalent, way of obtaining Eq. (27) has been found by Stodolsky. Stodolsky's formula is given by 
[17, 18] 

P = V X P - L>Pt, (37) 

where P is the "polarization vector", V is the "magnetic field", D is the damping parameter determined by the 
scattering amplitudes on the background {i.e., the environment), and Pt is the "transverse" part of P. In a two-level 
system, the length of P measures the degree of coherence: |P| = 1 corresponds to a pure state, < |P| < 1 to 
some degree of incoherence, and |P| = to the completely mixed or incoherent state [17], i.e., the loss of coherence 
is given by the shrinking of |P|. The time evolution described by Eq. (37) is a precession around V, combined 
with the shrinking of jP| to zero. The final state corresponds to p = ^, where both neutrino flavor states, and 
Ufi, are equally populated, and with vanishing coherence between them. In our notation, P = p, V = k, and 
DpT = {J22=i |da|^) P-Z^"=i(da-p)da = rf^ (p - (p • k)k). Since here k = -e^, we have that DPt = (P{p-p3ez) = 

(fi{Pl,P2,0). 

Thus, we can now construct the time-dependent neutrino density matrix as 

p^t) = l{I, + p{t).a) = U S'^ " ^ ) ■ (38) 

y\ J 2^''^^^' ' 2\Pi{t)+iP2{t) 1-P3(i) J ^ ' 
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As usual, the neutrino flavor states l^a), where a = e,/i, are superpositions of the neutrino mass eigenstates \i^a), 
where a = 1, 2, i.e., 

2 2 

Wa) = 51 ^aaWa) = XI ^aaW"')' " = e, /i. (39) 

a— 1 a— 1 

If we represent the mass eigenstates by the vectors = (1,0) and \i'2) = (0,1), then the flavor states become 
|fe) = (cos6', sin6') and ji^^) = (— sin^, cos0). Furthermore, if we assume that the initial state of a neutrino is |i^e), 
i.e., the system is prepared in a pure Vg state, then the initial condition for the neutrino density matrix is 

p(0)^Pee^M(^'e|=(^3i^^,^3^ COS^ ^ J = - (/2 + p(0) • a) , (40) 

and it follows, using Eq. (38), that pi(0) = sin 2^, P2(0) = 0, and P3(0) = cos 2^. Similarly, if the initial state is 
then we have 

I > , I / sin^ — sin 6* cos 9\ ^ , 
^''''^l^'')<^''l=^-sinecos^ cosH j=l-^-. (41) 

Thus, the neutrino transition probabilities with decoherence effects read 

Pee(<) = tr{\ue){ve\p{t)) = ]^[l + sin29pi{t)+cos2ep:i{t)], (42) 
Pe^(t) = tr {\Vf,){vM'^)) = ^ [1 - sin2^pi(i) - cos2eps{t)] . (43) 
Moreover, as in the neutrino oscillation formalism, we have the relation between the transition probabilities 

Pe^(i) = 1 - Pee(t) = 1 - P^^(t) = f^e(i). (44) 

since 

1 = trp(f) = tr {lp{t)) = tr ( ^ W){i^a\p{t) J = tr {\vg){vg\p{t)) + tr {\v ^,) {i^ M^)) = Pee{t) + Pe^,{t). (45) 

Inserting Eqs. (34) - (36) with the above initial conditions into Eq. (43) gives 

Pei^ {t) = ^ sin^ 26* (l - e-'^'* cos ktj . (46) 

Furthermore, inserting t ~ L and k = Am^/(2£'), we thus obtain the neutrino decoherence formula for two neutrino 
flavors as 

Pe.(i, E) = \ sin^ 26 (l - e-''^ cos ^) . (47) 
When Am^ = (no neutrino oscillations), we find the pure neutrino decoherence formula 

Pe^(L) = isin2 2^(l-e-'^'^). (48) 

Note that this formula is explicitly independent of the neutrino energy E. Other interesting limits are obtained when 
d = Q (no neutrino decoherence) 

Pe^(i,E) = -sin2 2^(^l-cos^^j , (49) 
i.e., the well-known Pontecorvo neutrino oscillation formula (c/. Eq. (8)) and when d ^ oo (or Am^ » 2E/L) 

P,^ = ^siv?2e, (50) 

i.e., the classical transition probability formula (c/. Eq. (7)), which is independent of both the neutrino path length 
L and the neutrino energy E. The effects of both neutrino oscillations and neutrino decoherence are said to be 
completely washed out. 
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III. EQUIVALENCE BETWEEN THE TWO MODELS 

In the previous section, we have seen that the two formaUsms, neutrino oscillations and neutrino decoherence, have 
the same dependence for the transition probabilities. In both scenarios, the oscillation factor cos is damped by 

an exponential factor (e^^"^ (^"^ ) or e^"^^^). 

Comparing Eqs. (13) and (47) with each other, we obtain 

2a^Am^)^ =cFL, (51) 
i.e., the decoherence parameter d is related to the damping parameter a as 



V2Am^ 

= ^ (52) 

which means that the two investigated models must be equivalent if this condition is fulfilled. The units of the 
decoherence parameter d is [d] = cV^^^. This investigation could of course be extended to n neutrino flavors, but we 
will not do this here. 

Next, we will now try to estimate the decoherence term (or the decoherence parameter d). The damping 
parameter a can be written as 



a = -^r, where r = ^ + ^ (pessimistic) or r = (^)^ + (^)^ (optimistic). (53) 



L 
4E 

Inserting Eq. (53) into Eq. (52) yields 



or 



Thus, we obtain the following estimate for the decoherence term 

Vocd^^ (56) 

if we assume that = 0(1). 

Recently, other estimates of the decoherence term 2?[p] have been found by Lisi et al. [3] and Adler [6]. Their 
different estimates are 

V ^ (57) 

Afpianck -^Planck 

and 

(58) 

Mpianck E Mpianck 

respectively, where again V ex d"^, — E, AHm — = k, and Mpianck is the Planck mass scale. They 

argue that the decoherence could be due to e.g. quantum gravity. We, on the other hand, have here argued more 
phenomenologically. 

In the next section, we are going to estimate the decoherence parameter for atmospheric neutrinos. 

IV. ESTIMATION OF THE DECOHERENCE PARAMETER FOR ATMOSPHERIC NEUTRINOS 
The path length for atmospheric neutrinos, which is easily obtained from geometrical considerations, is given by 

L = L(cos 1?) = cos2 7? + 2i?erf + - ij^ cos (59) 
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where d is the zenith angle, i?@ is the radius of the Earth (i?® ~ 6400 km), and d is the typical altitude of the 
production point of atmospheric neutrinos above the surface of the Earth {d ~ 10 km). The uncertainty in the path 
length is mainly determined by A cos ^? and therefore 



AL 



5L(cost?) 



9cosi? 



D 7" A r D 

Acosi?= , ® .Acosi?, i.e., =^ = , ® =Acosi9. (60) 

L + i?eCos^ L ^Rlcos^^ + 2R^d + d^ 

The uncertainty in the neutrino energy for the Super-Kamiokande experiment [19] as well as for the future MONOLITH 
experiment [20] is of the order of magnitude 

AE^E, i.e., ^ = 0{1), (61) 

even though the energy resolution is better for the MONOLITH than Super-Kamiokande they are of the same order 
of magnitude. 

Thus, using Eqs. (5), (60), and (61) as well as inserting numerical values, we have the following upper bound for 
the damping parameter (Tatm. for up-going (cos^? = —0.95 and Acos^? = 0.1 => L ~ 12000km and AL/L ~ 0.11) 
atmospheric neutrinos 

o-atm. ^ 3.0 • lO"'^ m/eV ~ 6.0 • 10"" m^ ~ 1.5 • lO'"' eV'^ (62) 

From two flavor neutrino oscillations analyses of Super-Kamiokande data, assuming lyn-iyr oscillations, the atmo- 
spheric mass squared difference Am^^j^ has been measured to be Am^^j^ ~ 3.2 • 10~^ eV^ (and the mixing angle 
^atm. — 45°) [19]. This means that the decoherence parameter for atmospheric neutrinos is 



rfatxn. = ^^"'■atn.^ < 2.8 • lO'^ eV^^ (63) 

VL(cos^ = -0.95) 
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or 

dltm. ^ 7.9 • 10~i^ eV = 7.9 • lO'^^ GeV - IQ-^" GeV. (64) 

Recently, Lisi et al. have discussed three scenarios for the decoherence parameter on the form d^^^ = 70 (cfv)"' 
where 70 is a constant and n = —1,0,2 [3]. Our result is in general agreement with the ones obtained by Lisi et al. 
{i.e., they are all comparable within some orders of magnitude), which are d^tm < 0.9-10~^^GeV at 90% C.L. (n ~ 2), 
d^^^ < 3.5 • 10-23 GeV at 90% C.L. [n = 0), dl^^ < 4.1 • lO^^s GeV at 95% C.L. (n = 0), dl^^ < 5.5 • 10-^3 GeV at 
99%'C.L. (n = 0), and d^^^ < 2 • lO^^i GeV at 90% (n = -1) [3]. They found that d^^^ oc E-^ is favored, which 
is natural from our point of view, since this would mean that the dependence of the exponentially damping factor in 
the transition probabilities is on the form e~"^^^ , where a is some constant, i.e., it is exponentially decreasing with 
the L/ E dependence. They also fitted the data with a pure neutrino decoherence model (Ato^^j^^ = 0) and obtained 
<^atm. = 1-2 • 10"^^ GeV (n = —1) [3]. In all their fits, they obtained maximal mixing, i.e., ^atm. = 45°, as their best 
fit value. 



V. SUMMARY AND CONCLUSIONS 

In conclusion, we have shown, in this paper, that a Gaussian averaged neutrino oscillation model and a neutrino 
decoherence model are equivalent if Eq. (52) is satisfied, i.e., d = Our upper bound on rfatm. ~ lO-^^GeV 

then immediately follows from the results of the Super-Kamiokande analyses, i.e., the neutrino decoherence model is 
just a reparameterization of the Gaussian averaged neutrino oscillation model or vice versa. We have also estimated 

the decoherence term to be 2? oc li^ ~ ^ , which is different from earlier results including the Planck mass scale 

Afpianck and motivations that decoherence could be induced by new physics beyond the Standard Model such as e.g. 
quantum gravity. 
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